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Abstract 

The far-infrared absorption of triangular and square two-dimensional quan- 
tum dots is studied by means of time simulations of the density oscillations 
within the time-dependent local-spin-density approximation. The absorption 
is spatially analyzed using a local-response theory that allows the identifica- 
tion of corner and side modes in the geometric nanostructures. The evolution 
with a vertical magnetic field of varying intensity is also discussed. 
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I. INTRODUCTION 



The quantum-dot field is recently focussing much attention on the analysis of symmetry 
unrestricted nanostructures, both from the experimental and theoretical point of view. An 
important issue is the influence of the dot geometry on the edge modes that manifest in 
the far-infrared absorption. It is well known that parabolically confined dots can only 
absorb energy at the frequencies of the confining parabolas (generalized Kohn's theorem 
However, either departures from a pure radial quadratic law or the introduction of 
an additional angular dependence in the confining potential may cause the fragmentation of 
the Kohn modes and the emergence of a more complicated pattern. 

At first, the great majority of theoretical approaches to quantum-dot absorption were 
restricted to the circular symmetry case, solving the random-phase- approximation (RPA) 
equations by fully exploiting the analytical angular dependence of the single-particle or- 
bitals Although very successful, this scheme is not easily applicable to symmetry- 

unrestricted nanostructures and different approaches have been proposed. In Ref. 
Magnusdottir and Gudmundsson address the problem within the Hartree theory by solv- 
ing the resulting RPA matrices in the single-particle basis of the circular problem. The 
enormous matrix dimensions make the calculation feasible only for a rather small electron 
number. Alternatively, Ullrich and Vignale [0 have used a perturbative method to describe 
the absorption of elliptic dots. A different approach was used by two of us in Ref. |^, con- 
sisting in the time simulation of the mean field density oscillations that follow an appropriate 
initial perturbation and from which one may extract the oscillation frequencies. This real 
time calculations do not rely on any symmetry restriction having also been used to analyze 
excitations in deformed atomic nuclei and metal clusters |[9|,p!0[]. 

In this work we concentrate on the description of the absorption in triangular and 
square shaped quantum dots using the time-dependent local-spin-density approximation 
(TDLSDA). The extension to low magnetic fields within this theory is considered as in Ref. 
1^. We show that the deformation leads to a departure from the Kohn modes which is 
peculiar to the shape. At low magnetic fields the spectra possess a rich structure being 
possible to identify the presence of peak anticrossings associated with the deviation from 



circular shape, similar to those found experimentally by Demel et al. To further clarify 



the origin of the fragmentation, an analysis of the local absorption at different points of the 
nanostructure is given. This allows to associate each excitation with a particular absorption 
pattern, thus giving a more physical interpretation of the peaks. In particular, a distinction 
between corner and side modes can be established in this way. The magnetic field is shown 
to introduce an angular average of the patterns leading, for increasing intensity, to sizeable 
shape modifications of the absorption distribution. 

Section II of the paper is devoted to the characterization of the ground state for the 
square and triangular quantum dots. In Sec. Ill we address the description of the far 
infrared absorption presenting in different subsections the formalism and global spectrum 
(III. A), the local absorption (III.B) and the angular average and magnetic evolution of the 
absorption patterns (III.C). The conclusions are given in Sec. IV. 
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II. SYSTEM CHARACTERIZATION 



We focus our attention on semiconductor nanostructures of similar sizes and different 
polygonal shapes, in order to discuss the differences introduced by the geometry. It is 
assumed that the electron positions are restricted to the xy plane where, in addition, an 
external potential is confining them to an electronic island. Different shapes are obtained 
by using the following potential, proposed by Magniisdottir and Gudmundsson 0], 

V;.i(r, 9) = ]^ujy{l + apCo^ipO)) , (1) 

where we have used planar polar coordinates (r, 9) and the standard effective-atomic-unit 
system [|12|. The potential parameter o^o gives the parabola curvature while p and deter- 
mine the deviation from circular symmetry. In this work we shall restrict to triangular and 
square shapes. These are the two simplest geometries that deviate from a pure parabolic 
behaviour and whose absorption is therefore not fixed by the generalized Kohn's theorem. A 
common parabola frequency ujq = 0.5 a.u. has been chosen both for the triangle and square 
while the deformation parameters have been fixed to as = 0.31 and = 0.15 respectively, 
in order to obtain a proper triangle and box shape, as can be seen in Fig. 1. 

We have described the electronic structure within density-functional theory in the LSDA. 
This approach has been used to address a variety of nonuniform electronic systems — such 
as atoms, molecules and clusters — and quantum dots are not an exception for many authors 
have used this scheme to analyze dot properties [p!^ -[T6[| . Obviously, for our present purposes 



a symmetry unrestricted solver is required. This can be obtained with different techniques 
|7| J15|JT6[| such as diagonalization in a plane-wave basis or perturbative treatements based 



on a circular solution. We have resorted to the method developed by two of us in Ref . p , 
based on the uniform discretization of the xy plane and the associated transformation of the 
diferential Kohn-Sham equations to be solved with iterative techniques. This method has 
also been applied by us to the study of quantum rings and dot molecules 



In what follows we have restricted the present study to nanostructures containing an in- 
termediate number of electrons (A^ = 12) for which LSDA is expected to provide a reasonable 
description of electron-electron interactions. Figure 1 shows that the electronic distribution 
closely follows the shape of the external potential, clearly of triangular or square form. How- 
ever, there are also conspicuous density peaks localized either at the corners or the interior 
of the structure. They are of course attributed to the formation of delocalized electronic 
orbitals, and to quantum effects associated with the presence of boundaries in the potential, 
similar to the well-known Friedel oscillations at a metal surface. 



III. FIR ABSORPTION 

A. Formalism and spectra 

The description of excitations in deformed nanostructures constitutes a highly nontrivial 
task because of the lack of symmetry which does not allow the analytic integration of angular 
variables as in circular systems. In the present work we have used the TDLSDA to obtain 
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the time evolution, following an initial perturbation, of relevant expectation values. The 
single-particle orbitals {(fi{r)} evolve in time as 



d 

i-^(fiinir, t) = h^[p, m] (/?i^(r, t) , (2) 

where 1] =t, i is the spin index, and total density and magnetization are given in terms 
of the spin densities p-^ir) = J2i lv^jr?(r)P, P = Pi + Pi and m = p| — pi, respectively. 
The Hamiltonian hrj in equation (0) contains, besides of kinetic energy, the confining and 
yi'^°^f) (^y) , Hartree t>^^^(r) = / dr' p{r')/\r — r'\ exchange-correlation ^^^'^^(r) = -^8xc{p,rn) 
potentials. The exchange-correlation energy density S^dp, has been described as in Refs. 
I^J^. The Zeeman energy ±{l/2)g*pBB [0] is also included in Eq. ||. 

An initial perturbation of the ground state orbitals v^'(r) = V(p{r) models the inter- 
action with the external field. For dipole excitations, corresponding to FIR absorption in 
nanostructures, the unitary operator V is represented by a rigid translation in an arbitrary 
direction e, specifically 

V = exp [iX e ■ p] , (3) 

where the A parameter has been taken small enough in order to keep the system response in 
the linear regime. After the initial excitation we keep track of (r)(t) from where a frequency 
analysis provides the absorption energies and their associated intensities characterized 
by the strength function Se{uj), 

PeH =e- J dte"'' (r)(t) , 

S'eH = |PeM|. (4) 

Usually, one is interested in an average over different oscillation directions, corresponding to 
the absorption of nonpolarized radiation. This implies to average Eq. (4) over all e directions. 
Taking into account that in the small amplitude limit the strength in a particular direction 
only depends on the initial amplitude given in that direction the averaged strength may be 
obtained as pO 
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[ " de |cos(^)2p^(cu) + sm{ef Vy{uj) I . (5) 

^ 
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It can easily be shown that for geometries with two or more equivalent symmetry axes S^yuj) 
actually does not depend on e and therefore it coincides with Sav{^^)- 

Figure 2 displays the averaged spectrum obtained at i? = for the square and triangle. 
Notice that both cases show a two peak spectrum within LSDA, although the interpeak 
separation is much higher for the triangle than for the square. This can be attributed to 
the higher level separation in the triangle. In fact, in the triangular case the single-particle 
levels exhibit a shell-like structure with gaps between groups of close levels, while for the 
square level straggling is somewhat higher. This reflects in the distribution of particle-hole 
transitions also shown in Fig. 2. 
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The evolution of the absorption spectra with the intensity of the static magnetic field is 
shown in Fig. 3. In general, the magnetic field induces a splitting of the B = peaks, as is 
well known for the circular parabola [Q] where the two frequencies are u!± = ^Ju^ + u'^/A ± 
a;c/2. For the square and triangular dots the situation is more involved, since the 5-induced 
splitting goes along with geometry effects. Figure 3 indicates that there is a repulsion 
amongst the branches that originate from the B = modes, showing a zone of avoided 
crossings or anticrossings of the collective peaks. These anticrossing areas are very sensitive 
to the geometry, showing complicated patterns with more than just two peaks. Indeed, for 
the square it roughly covers the interval 0.3 — 0.7 T while for the triangle this region is 
more extended lying between 0.5 — 2 T. A similar difference is found when comparing the 
magnetic field Bpoi at which the ground state spin polarizes from 5* = to 5 = 1. The 
square first transition lies at Bpoi = 0.5 T while for the triangle it is Bpoi = 2.7 T. We also 
notice that the minimum energy gap between the two main branches is much smaller for the 
square than for the triangle. 



B. The local absorption 

To further analyze the origin of the fragmented absorption in squares and triangles we 
have computed the local absorption in a way similar to that used in Ref. |^ in the context of 



excitonic states in semiconductors. This technique allows us to study the global absorption 
spectrum taking into account the different contributions from each spatial point, depending 
on how the electronic density is oscillating in time in the vicinity of that point. In practice, we 
simulate this position-dependent probe by introducing a local weighting function (a narrow 
Gaussian centered on each point) and look for the oscillation frequencies of the convoluted 
density as a function of the probe position. After an initial dipole perturbation each local 
signal will evolve differently in time, thus manifesting the differences in absorption on a local 
scale. Figure 4 shows the local absorption for the two modes at i? = of the square and 
triangle when the system is excited along the x-axis (e = x). As a reference we also show 
the circular parabolic case, namely, ctp = in Eq. (0). 

The density in a circular parabola oscillates rigidly at the parabola frequency (Kohn's 
theorem) and from Fig. 4 we notice that the absorption indeed localizes on the edges along 
the oscillation direction (x-axis), as was to be expected. Turning now to the square geometry 
we notice that the dominanting mode at uj = 0.49, very close to the parabola frequency 
uo = 0.5, spatially locahzes along the edges in the polarization direction in a way quite 
similar to the circular case. On the contrary, the lower square mode at = 0.42 very clearly 
deviates from this situation showing four maxima at the corners of the nanostructure. This 
result suggests a classification of the two excitations as side and corner modes, respectively. 
This classification is confirmed in the case of the triangle, where we see that the u = 0.73 
mode basically absorbs on one of the triangle sides while the u = 0.44 absorbs on the 
opposite corner. The distinction between side and corner modes will be further commented 
in the next subsection where we shall present the angular average of the absorption patterns. 

A remarkable property seen in Fig. 4 regards the spatial self-organization of the absorp- 
tion modes. Indeed, we realize that both in the square and the triangle the absorption 
maxima for a given mode lie at regions where the absorption for the other mode is strongly 
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quenched. Wc may interpret this as a mode competition for an efficient distribution over the 
available nanostructure surface. 

C. Magnetic field dependence 

The magnetic field introduces a Coriolis-like force that makes the oscillation axis rotate in 
time. As a consequence, the initial chosen direction becomes unimportant since the magnetic 
field will automatically change it in time. We can take advantage of this effect to obtain 
the averaged absorption patterns of Fig. 4 over polarization directions by switching on a 
very small magnetic field. Figure 5 shows the angularly averaged patterns corresponding to 
the same modes displayed in Fig. 4. The distinction between corner and side modes is very 
clear in the case of the triangle: for u; — 0.44 the mode is of corner type and for u; — 0.73 
it is of side character. For the square this distinction is not as clear because the corners 
always absorb more than the side centers. Ncvertheles, the corner character is somehow 
higher for the low energy mode uj = 0.42 than for the high uj = 0.49 one; in fact the former 
shows vanishing absorption at the side centers while the latter has a much smoother angular 
dependence. 

When the magnetic field is increased additional fragmentation, with a complicated distri- 
bution of the absorption, appear both in the square and the triangle (cf. Fig. 3). Because of 
the very large number of possible cases, we shall only select the more illustrative examples in 
the anticrossing areas. Our aim is to discern the absorption pattern evolution with applied 
magnetic field and in particular, to check whether there is a relation between the energy 
repulsion of modes and their spatial localization. Figures 6 and 7 show the evolution with B 
of the pattern for the two branches with the higher absorption, that repel each other, in the 
square and triangular dot, respectively. We see in Fig. 6 that indeed for the square a clear 
localization takes place a.t B = 0.4 T (the center of the anticrossing region) at the corners 
for low energy branch and at the sides for the high energy one. In the triangle, a more 
conspicuous side and corner character of the branches is retained for all magnetic fields due, 
as discussed above, to the existence of a greater energy gap between the different modes. 
Figure 7 shows the patterns for = 1 T as well as the corresponding energy in the spectral 
range. The three lower patterns give rise, when added according to their relative strength, 
to the global strength distribution displayed in the rightmost panel which exhibits a clear 
corner identity. Conversely, the upper branch possess a markedly side character. The same 
same trends have been obtained in the full range of magnetic fields corresponding to Fig. 3 
for the triangle. 

IV. CONCLUSIONS 

The FIR absorption of a quantum dot with a triangular or square shape has been ana- 
lyzed taking into account the effect of an applied static magnetic field, by means of real-time 
simulations in a spatial mesh within the LSDA formalism. The fragmentation of the ab- 
sorption due to the noncircular shape has been emphasized. In particular, the existence of a 
rich absorption distribution for intermediate magnetic fields, associated with the repulsion 
between modes in the so-called anticrossing energy interval has been analyzed. 
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The absorption pattern corresponding to the several excitation modes, for a given mag- 
netic field, has been obtained from the local density oscillations. These patterns allow us to 
easily interpret the physical character of each mode. In particular, we have identified the 
existence of corner and side modes in the nanostructure. In the anticrossing intervals the 
absorption patterns of different modes do not significantly overlap, a result that we atribute 
to a mode competition for the available dot surface as well as a mode spatial repulsion, 
consistent with the energy repulsion leading to the anticrossing effect. 
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FIGURES 



FIG. 1. Density distributions for the square and triangular dots with 12 electrons considered in 
this work (see Sec. II). Upper plots display the densities as three dimensional surfaces while lower 
plots show the corresponding contour lines. 

FIG. 2. Absorption spectra for the square and triangle. Solid line corresponds to the TDLSDA 
while the dashed line shows the non-interacting response. For clarity the non-interacting response 
has been slightly shifted in the vertical direction. The right plots show the relative position of the 
lowest energy levels in each case with the arrow indicating the position of the highest-occuppied 
orbital. 

FIG. 3. Magnetic field evolution of the absorption spectrum of the square and triangle shown 
in a grey colour scale, with black and white indicating high and low aborption levels, respectively. 

FIG. 4. Local absorption patterns in the circle, square and triangle for an oscillation along the 
X-axis and the energies indicated in each panel. 

FIG. 5. Average over polarization directions of the same modes of Fig. 4 obtained by switching 
on an small vertical magnetic field. 

FIG. 6. Magnetic field evolution of selected absorption patterns within the anticrossing region 
of the square dot. The magnetic field and energy allow to exactly identify what mode of Fig. 3 is 
being considered in each case. 

FIG. 7. Spatial absorption patterns corresponding to different modes of the triangular nanos- 
tructure for S = 1 T. The arrows join each spatial pattern to its corresponding excitation energy 
in the spectral interval. See text for more details. 
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